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Abstract. In recent years, many difficulties appeared when taking into account the inherent 
stochastic behavior of neurons and voltage-dependent ion channels in Hodgking-Huxley type 
models. In particular, an open problem for a stochastic model of cerebellar granule cell 
excitability was to ensure that the values of the gating variables remain within the unit 
interval. In this paper, we provide an answer to this modeling issue and obtain a class of 
viable stochastic models. We select the stochastic models thanks to a general criterion for 
the flow invariance of rectangular subsets under systems of stochastic differential equations. 
We formulate explicit necessary and sufficient conditions, that are valid for both, Ito's and 
Stratonovich's interpretation of stochastic differential equations, improving a previous result 
obtained by A. Milian [A.Milian, Coll. Math. 1995] in the Ito case. These invariance criteria 
allow to validate stochastic models in many applications. To illustrate our results we present 
numerical simulations for a stochastic Hodgkin-Huxley model. 



1. Introduction 

In recent years, a great deal of activities has been devoted to develop models of neuronal 
excitability that take into account the intrinsic stochastic bioelectrical activity of neurons 
(see [11], [3]). In [11], the authors apply Ito's theory of stochastic differential equations and 
propose a stochastic model which reproduces the irregular electrophysiological activity of an 
in vitro granule cell (see [llj and Figure 6 p. 7). A particular case of this model is a stochastic 
version of the classical Hodgkin-Huxley model (see [I], [3]). However, as already indicated by 
the authors, the model suffers severe difficulties ([11] p. 4 and p. 10): Undesired values were 
observed for the gating variables that are supposed to take values within the unit interval. 
The solution of this modeling difficulty is mentioned as a challenge for future work ( jllj p. 10). 
Similar problems also occurred for the stochastic Hodgkin-Huxley model in [3] (see p. 2071). 
In this article, we provide an answer to this problem and obtain a family of viable stochastic 
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models for cerebellar granule cell excitability. The admissible models are derived from a gen- 
eral invariance theorem for systems of stochastic differential equations. 

We formulate invariance results in a general setting that allow to validate stochastic models 
in many applications. When the solutions of a given system of stochastic differential equa- 
tions describe quantities that necessarily take values within a certain range, the problem can 
be mathematically analyzed by studying the flow invariance of rectangular subsets of the 
euclidean space. To be more precise, we consider systems of Ito differential equations of the 
form 

dX(t) = f{t, X(t))dt + g(t, X(t))dW{t), 

where the process X is vector- valued. For the concrete formulation we refer to Section [2] We 
characterize the class of functions / and g that lead to viable stochastic models and formulate 
explicit necessary and sufficient conditions that can be directly checked and easily verified 
in applications. We further show that the invariance theorems are valid for both, Ito's and 
Stratonovich's interpretation of stochastic differential equations. The result for the Ito case 
was previously obtained by A. Milian in [7J. As discussed in [12] , in a concrete application 
it is generally not easy to decide which interpretation should be applied. Our results show 
that the qualitative behavior of solutions regarding non-negativity and boundedness is inde- 
pendent of Ito's or Stratonovich's interpretation. Other properties of the solutions, however, 
may strongly depend on the choice of the interpretation (see [8]). 

The outline of our paper is as follows: In Section [2] we introduce the class of stochastic sys- 
tems we study and formulate general invariance criteria for systems of stochastic differential 
equations. We then apply the results to obtain viable stochastic models for cerebellar granule 
cell excitability in Section [3} In Section [4] we present numerical simulations to illustrate the 
model behaviour. Finally, in Section [5] we recall the results obtained by A. Milian in [7] and 
present the proofs of the invariance theorems. 



2. Invariance Criteria for Stochastic Differential Equations 

Let (Q, J-, P) be a probability space with a right-continuous increasing family F = (J-t)t>o of 
sub-cr-fields of T each containing all sets of P-measure zero. We consider systems of stochastic 
Ito equations of the form 

X(t) = X + f f(s, X(s))ds + ! g(s, X(s))dW(s), t G [t , oo[, (1) 

J to J to 

where / = [/•] : [0,oo[xM m R m is Borel-measurable, and g = [gij] : [0,oo[xM m R mxr 
is a Borel-measurable mapping into the set of all M mxr -matrices, i = 1, . . . ,m, j = 1, . . . , r. 
Furthermore, W : [0,oo[xf2 — > W denotes an r-dimensional F- Wiener process, the initial 
time to is non-negative and X$ £ M. m is the given initial data. 

The stochastic integral equations ([!]) are commonly written as system of stochastic Ito differ- 
ential equations, 



dX(t) = f(t,X(t))dt + g(t,X(t))dW(t), t€[t ,oo[, 
X(to) = Xq, 
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where the function g represents the stochastic perturbation and / the deterministic part. 
Indeed, if g = we obtain the corresponding unperturbed deterministic system of ODEs. 

In the sequel, we denote by (f,g) stochastic initial value problems of the form ([!]). We aim 
at formulating explicit necessary and sufficient conditions on the functions / and g for the 
non-negativity and boundedness of solutions. The conditions can directly be verified and 
allow to explicitly characterize the class of admissible models in applications. In Section [5] 
we deduce our main theorems from a more general result about the stochastic invariance of 
polyhedral subsets of M m . However, in applications the non- negativity and boundedness of 
solutions are the most relevant modeling issues. 

Since our aim is not to establish the well-posedness of the stochastic initial value problem but 
to study the qualitative behavior of solutions, we assume that for every initial time to > 
and initial data Xq G IR m there exists a unique solution of the stochastic problem Q. 

Definition 1. We say that the subset K C K m is invariant for the stochastic system ( f,g) 
if for every initial data Xq G K and initial time to > the corresponding solution X(t), 
t > to, satisfies 

P{{X(t)EK, i€[to,oo[}) = l. 

The following theorem characterizes the class of functions / and g such that the stochastic 
system (/, g) preserves the positivity of solutions. That is, solutions emanating from non- 
negative initial data (almost surely) remain non- negative as long as they exist. 

Theorem 2. Let I C {1, . . . , m} be a non-empty subset. Then, the set 

K + := {x = (xi, . . . , x m ) G M. m , x { > 0, i G /} 

is invariant for the stochastic system (/, g) if and only if 

fi(t,%) > for x G K + such that x,- t = 0, 
gij(t, x) = for x G K + such that = 0, j = 1, . . . , r, 

for all t > and i G I. 

This result applies independent of ltd 's or Stratonovich 's interpretation of stochastic differ- 
ential equations. 

The solutions of mathematical models often describe quantities that necessarily take values 
within a particular range. We next formulate a criterion for the invariance of rectangular 
subsets of the phase space W 71 . 

Theorem 3. Let I C {l,...,m} be a non-empty subset and ai,b{ G M such that bi > a{. 
Then, the set 

K := {x G R m : a { < x { < b h i G /} 
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is invariant for the stochastic system (/, g) if and only if 

fi(t,x) > for x £ K such that X{ = a., 
fi(t,x) < for x £ K such that Xi = bi, 
gij(t, x) = for x G K such that xi £ {aj, bi}, j = 1, . . . , r, 

for all t > and i £ I. 

This result is valid independent of ltd 's or Stratonovich 's interpretation. 

If we apply Theorem [2] or Theorem [3] to the corresponding unperturbed deterministic system 
(/, 0) we recover the well-known tangential condition for systems of ODEs, which is necessary 
and sufficient for the flow invariance of subsets of IR m (see |13j or [S]). 

Finally, we formulate a criterion for the validity of comparison principles for the solutions of 
stochastic systems. A. Milian stated the following theorem for systems of ltd equations in [7J. 
We recall her result and show that it remains valid if we apply Stratonovich's interpretation 
of stochastic differential equations. 

Theorem 4. Let I be a non-empty subset of {1, . . . ,m}. We assume that (f,g) and (f,g) 
are stochastic systems of the form with given initial data Xq,Yq G M. m , and denote by X 
and Y the corresponding solutions. Then, the following statements are equivalent: 

(a) For all tQ > and i £ I, if the initial data satisfy (Xq)^ > (lo)i> then 

P{{X t {t) > Y t (t), t £ [to, oo[, i G /}) = 1. 

(b) The functions f,f,g and g satisfy 

fi(t,x) > fi(t,y), t>0, 
9ij(t,x) = gij(t,y), t > 0, j = l,...,r, 

for all i £ I and x,y £ W 11 such that Xi = yi and x^ > for k £ I. 

3. The Validation of a Stochastic Hodgkin-Huxley Type Model 

A stochastic model for cerebellar granule cell excitability was proposed and numerically stud- 
ied in As many biophysical models of neurons it is based on the well-known deterministic 
Hodgkin-Huxley formalism [1], which qualitatively describes the conduction and excitation 
in nerves. Such models are commonly formulated as systems of deterministic ODEs. The 
behavior of neurons and voltage-dependent ion channels, however, is known to be stochastic 
in nature, which motivates the modeling approach in [H]. The mathematical model is formu- 
lated as system of stochastic differential equations for the dependent model variables Xj, which 
represent the gating variables for the specific ion channels, the transmembrane potential V 
and the intracellular calcium concentration C, 

dxi = fi(V,Xi)dt + aidW.it), i = l,...,8, (2) 

dx 9 = f 9 (V, C, x 9 )dt + cr 9 dW 9 (t), 
dV = F(t,V, xi, . . . , x 9 )dt, 
dC = G(V,C, x 7 ,x 8 )dt, 
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and the rate functions for activation on and inactivation /3j are continuous and positive. The 
stochastic differential equations are interpreted in the sense of Ito, Wi(t),t > 0, denote stan- 
dard scalar Wiener processes, dW{ the corresponding Ito differentials, and the parameters <jj 
are positive and constant, % = 1, ... ,9. For the concrete form of the interaction functions F 
and G and the complete description of the model we refer to 

This model extends a previous deterministic model for cerebellar granule cell excitability by 
adding the stochastic terms (TjdWi(i) in the governing equations for the gating variables Xi, 
i = 1, ... ,9. Ion channel stochasticity has been detected experimentally and is due to the 
thermal interaction of molecules constituting an ion channel. It can be observed as random 
opening and closing of an ion channel at an experimentally fixed membrane potential (see 
PI] and also [3]). 

The gating variables Xi describe the opening and closing rates of the specific ion channels 
and necessarily take values within the interval [0, 1]. While the corresponding unperturbed 
deterministic model, where G{ = 0, i = 1, . . . , 9, certainly ensures this property, it cannot be 
guaranteed by the stochastic model Q: 

The parameters cij, i = 1,...,9, which take into account the intensity of the stochastic 
perturbations, were taken to be constant in the model and the simulations presented in 
The necessity to carefully choose these parameters was indicated. In particular, undesired 
values of the gating variables were observed and discussed, it was stressed that this modeling 
issue needed to be solved and highlighted as a challenge for future work (see [IT] , p. 4 and 
p. 10). Similar difficulties also occurred for the stochastic Hodgkin-Huxley model developed 
in [3] (see p. 2071). Our results show that independent of the choice of the parameters Oi the 
invariance of the unit interval cannot be guaranteed by the stochastic model Q if we take 
these parameters to be constant. Indeed, the conditions on the stochastic perturbations in 
Theorem K3| applied to the model Q and the invariant subset 



are never satisfied. 

We obtain viable stochastic models if we replace the constants Oi by appropriate functions 
gi, that ensure the desired invariance of the unit interval. To be more precise, we propose to 
consider models of the form 



K = {yeR n , 0<yi<l,i 



= 1 



-..,9} 



where x = (xi, . 



dxi 
dxg 
dV 



dC 



fi{V, Xi )dt + 9i {t, V, C, x) dW l (t), i = 1, . . . , 8, 

f 9 (V, C, x 9 )dt + gg(t, V, C, x) dW 9 (t), 
F(t,V,x)dt, 
G(V, C, X7, xg)dt, 

xg), and the stochastic perturbations gi : [OjOofxIR 11 — > M. satisfy 



(3) 



gi(t, y) = for y G K such that yi G {0, 1}, 
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for all t > and i = 1, . . . , 9. 



Proposition 5. The modified stochastic model ensures that the gating variables X{ take 
values within the interval [0,1], for alii = 1,...,9. This is valid for ltd 's and for Stratonovich's 
interpretation of the stochastic differential equations. 

Proof. The statement is a direct consequence of Theorem [3] since the interaction functions 
fi and stochastic perturbations gi in the governing equations for the gating variables Xi, 
i = 1, . . . , 9, in the modified stochastic model Q satisfy the required conditions. □ 

One possible choice for the stochastic perturbations are functions of the form 

gi(xi) = UiXi(l - Xi), 
with constants o"j E M + , i = 1, . . . ,9. 



4. Numerical Simulations 

To illustrate our results we present numerical simulations for a simplified version of the 
stochastic model discussed in the previous section and consider a stochastic version of the 
classical Hodgkin-Huxley model [4] Despite its simplicity, the deterministic Hodgkin-Huxley 
model has always been playing a very important role in the study of neuron excitability ([6]). 
However, stochasticity should be included in the model to take into account the stochastic 
behavior of the ion channel kinetics (see [6], p. 558 and p. 559). 




5 10 15 20 25 5 10 15 20 25 



(a) Gating Variables (b) Voltage 

Figure 1. Deterministic Model 



The model we consider is formulated as system of ordinary differential equations for the 
dependent model variables Xi, i = 1,2,3, that represent the gating variables for the specific 
ion channels, and the voltage V, 

ai (V)(l - Xi ) - Pi(V)xi, (4) 

^[1 - 9Nax\xz(y - E Na ) - g K x\(y - E K ) - g L (V - E L )]. 



dxi 

dt 
dV 

~dt 
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The rate functions for activation and inactivation are given by 
0.1(V + 35) 



at(V) 
a 2 (V) 



l-exp(-^P)' 
0.01(^ + 50) 



V+50 \ ' 
10 > 



1 — exp( 

a 3 (V) = 0.07exp(-0.05(F + 60)), 
and the parameter values by 



C = 0.01- 



cm/ 



I = 0.1 mV, 
(see @). 



mS 

9Na = 1-2 

Ejva = 55.17 ml/, 



01 (V) = 4.0exp(-0.0556(y + 60)), 
V + 60 N 



p 2 (V) = 0.125 exp(- 



80 



1 



l + exp(-0.1(V + 30))' 



5x = 0.36 — 

cm z 

E K = -72.14 mF, 



9L = 0.03 — 
E L = -49.42 ml/, 




tage Change [o.r Hodgkln-Hir*li 



(a) Gating Variables: a = 0.1 



1 V 



(c) Gating Variables: a = 0.5 




(b) Voltage: a = 0.1 



Voliage Change far Hadgkin-Hu: 




(d) Voltage: cr = 0.5 



Figure 2. Stochastic Model: Additive Noise 



We illustrate the model behavior in Figure [TJ The gating variables Xj, i = 1,2, 3, describe the 
opening and closing rates of the specific ion channels and necessarily take values within the 
interval [0, 1]. The deterministic model Q certainly ensures this property. 



Following the modeling approach in |llj . we may extend the deterministic model by adding 
the stochastic terms Cj dWi (t) in the governing equations for the gating variables Xi in the 
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model Q, which leads to the system of stochastic differential equations 

dxi{t) = ( ai (V(t))(l - Xi (t)) - Pi(V(t))xi(t))dt + <Ji dWi{t), i = 1,2, 3. (5) 

We interpret the stochastic differential equations in the sense of Ito, Wi(i),t > 0, denote stan- 
dard scalar Wiener processes, dWi the corresponding Ito differentials, and the parameters <7j 
are positive and constant, i = 1, . . . , 3. 




(a) Gating Variables: <r = 0.1 (b) Voltage: a = 0.1 




(c) Gating Variables: a = 0.5 (d) Voltage: a = 0.5 



Figure 3. Viable Stochastic Model: Ito's Interpretation 

Our results in Section [2] imply that the gating variables in the model ^ take undesired values 
outside of the unit interval. The simulations in Figure [2] illustrate this observation for differ- 
ent values of the parameter a = (Ti,i = 1,2,3. Here, we used the Euler-Maruyama method 
for the numerical implementation (see [5] and [H3]). We remark that Ito's and Stratonovich's 
interpretation yield the same solution for the stochastic model Q. 

We obtain viable stochastic models if we replace the constants cij in the equations for the 
gating variables by appropriate functions gi, i = 1,2,3, that ensure the desired invariance 
of the unit interval. To be more precise, we may consider stochastic models, where the 
determining equations for the gating variables x% are of the form 

dxi(t) = (a i (y(t))(l-x i (t))-ft(y(t))x i (0)dt + 5i(i^i(i)^2(t),a;3(i))dW" i (t), (6) 

for i = 1, 2, 3, and the stochastic perturbations gi : [0, oo[xM 3 — > M satisfy 

gi(t, xi, X2, £3) = for x\, X2, X3 6 [0, 1] such that X\ £ {0, 1}, 
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Voltage Change for Hadgkin-Hu: 











h 










V 





(a) Gating Variables: cr = 0.1 




(c) Gating Variables: tr = 0.5 




(b) Voltage: a = 0.1 




(d) Voltage: cr = 0.5 



Figure 4. Viable Stochastic Models: Stratonovich's Interpretation 
for all t > and i = 1, . . . , 3. 

Theorem [3] in Section [2] immediately implies that the gating variables in the modified sto- 
chastic model ([6]) take values within the interval [0,1], and that it is valid for Ito's and 
Stratonovich's interpretation of the stochastic differential equations. 



We illustrate the model behavior for stochastic perturbations of the form 

gi(t, x 1 ,x 2 , x 3 ) = axi(l - Xi) i = 1, . . . , 3, 

where the constant a > 0. The simulations in Figure [3] show the behavior of the solutions for 
Ito's interpretation of the stochastic system (pi), where we used the Euler-Maruyama method 
for the numerical implementation. Figure [^illustrates the model behavior when we apply 
Stratonovich's interpretation. In this case we applied the Euler-Heun method for the simula- 
tions (see [5] and [TO]). 



5. Proof of the Theorems 

Our proof is based on the main theorems obtained by A. Milian in [7]. We first recall her 
results, which are formulated for systems of stochastic ltd differential equations, and yield 
necessary and sufficient conditions for the stochastic viability of polyhedral subsets and the 
validity of comparison theorems. 
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Definition 6. A subset K C M m is said to possess the stochastic viability property with 
respect to the system (/, g) if for every initial data Xq G K and every t$ > § there exists a 
global solution of the initial value problem and the solution satisfies 

P({X(t)eK, te [to,oo[». 

For vectors a,n G M m we denote by 

H a , n ■= {x G R m , (x -a,n)> 0} 

the half-space determined by a and n, where (•, •) is the scalar product in M m . A polyhedron 
K in M m is a set of the form 

where / = {1, . . . , N} C N is a finite subset and a u , n u G M m , v G I. 
For the proof of the following results we refer to [7] . 

Theorem 7. Let K = f\e/ ffa^,n„ be a polyhedron in W n and suppose that the functions f 
and g satisfy the following conditions: 

(a) For every T > there exists a constant Ct > such that 

\\f{t,x)\\ 2 + \\g{t,x)\\ 2 < C T {1 + \\xf ) for all x G K,t G [0,T]. 

(b) For every T > i/iere exists a constant Ct > suc/i i/icrf 
||/(t,x)-/(t,y)|| + || 5 (t,x)-«7(t,y)|| <C T ||x-y|| /or a// x,yeK,te [0,T\. 

(c) For every x £ K the functions f(-,x) and g(-,x) are continuous on [0, oo[. 

Then, the set K possesses the stochastic viability property with respect to the system (/, g) if 
and only if for all u £ I and x G K such that {x — a u , n u ) = we have 

(f(x,t),n u ) >0, 

(9j(x,t),n„) = 0, j = l,...,r, 
for all t > 0, where gj is the j-th column of the matrix g = [gij]. 

Theorem 8. Let I be a non-empty subset of {1, . . . ,m} and suppose that for every T > 
there exists a constant Ct > such that 

(a) \\f{t,x)\\ 2 + \\g(t,x)\\ 2 <C T (l + \\x\\ 2 ) for all x G M m ,t G [0,T]. 

(b) + <C T ||x-y|| /or a// x,y £ W 71 , t G [0, T]. 

(c) For every x £ K the functions f(-,x) and g(-,x) are continuous on [0, oo[. 

We assume that the functions f and g satisfy the same conditions and denote the correspond- 
ing solutions of the stochastic systems (f,g) and (f,g~) by X and Y. Then, the following 
statements are equivalent: 

(i) For all t > 0, X = ((X )i, . . . , (X ) m ) G R m and Y = ((F )i, • ■ ■ , (*o)m) G W m 
such that (Xo)i > (Yo)i,i G /, the corresponding solutions satisfy 

P{{Xi(t)>Yi{t), i€l,t>to}) = l. 
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(ii) For all i £ I the functions f and g satisfy 

fi(t,x) > fi(t,y) fort>0, 
9ij(t, x) = g {j {t, y) for t > 0, j = 1, . . . , r, 

and all x = (x%, . . . ,x m ) G M m , y = (y±, ■ . . ,y m ) € M m such that Xk > yu, k 6 /, 
x% — y% ■ 

We remark that the assumptions (a)-(c) in Theorem [7] are imposed to guarantee the existence 
of solutions of the stochastic initial value problem (f,g). The more restrictive hypothesis 
in Theorem [8] imply the existence and uniqueness of solutions, and therefore, the stochastic 
viability of a subset is equivalent to the stochastic invariance with respect to the system (/, g). 

We will deduce our criteria from Milian's result and show that they are valid independent 
of Ito's and Statonovich's interpretation of stochastic differential equations. In the sequel, 
we use the symbol o dW(t) to indicate Stratonovich's interpretation. For convenience of the 
reader we recall the general conversion formula for systems of stochastic differential equations, 
which relates both interpretations (see [2], Section 6E): 

If we interpret the stochastic system ([!]) in the sense of Stratonovich, that is, X is the solution 
of the stochastic system 

dX(t) = f(t, X(t))dt + g(t, X(t)) o dW(t), 
then, X solves the system of Ito equations 



dX(t) 



f(t,X(t)) + h(t,X(t)) 



dt + g(t,X(t))dW(t), 



where the function h = [hi] : [0, oo[xM m — > W 71 is given by 



hi(t,x) = ^2^2~-^(t,x)g jk (t,x), i = l,...,m. (7) 



k=i j=i 3 



d Xj 

Proof of Theorem 2. Using polyhedral subsets of M m the positive cone can be represented as 

m 

K+ = f]H 0>ei , 

i=l 

where £ M. m denotes the origin and e^, i = 1, . . . , m, the standard orthonormal basis vectors 
in W" 1 . Since we a priori assume the existence and uniqueness of solutions of the stochastic 
initial value problem (f,g), the stochastic viability of the positive cone K + is equivalent 
to its invariance with respect to the system (f,g). Evaluating the necessary and sufficient 
conditions formulated in Theorem [7] immediately follows the result for Ito's interpretation 
of the stochastic system. We need to show that the statement remains valid if we apply 
Stratonovich's interpretation of stochastic differential equations. Let X be a solution of the 
Stratonovich equation 

dX(t) = f(t, X(t))dt + g(t, X) o dW(t). 

Then, the transformation formula implies that X solves the system of ltd equations (f,g) 
with modified interaction term / = / + ^h, where the function h is defined by the formula 
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0. We apply our previous result, which is valid for Ito's interpretation, to the stochastic 
system (/, g) and conclude that the positive cone is an invariant subset if and only if 

fi(t> x ) > x G K + such that x\ = 0, 

gij(t, x) = x G such that Xi = 0, j = 1, . . . , r, 

for all t > and i £ I,. The conditions on the stochastic perturbations yield the representation 
gi,j(t,x) = Xi J -^-(t,xi, . . . ,sxi, . . . ,x m )ds i G I,= 1, . . . ,r, (8) 

and it follows that the functions / and g satisfy the conditions in Theorem [2] if and only 
if the functions / and g fulfill these conditions. This observation concludes the proof for 
Stratonovich's interpretation. □ 

Proof of Theorem 3. We can represent the subset K in Theorem [3] as the finite intersection 
of polyhedral subsets 

K = f]{H ^nH e ^ ei ). 

Computing explicitly the necessary and sufficient conditions for the invariance of the subset 
K in Theorem [7] follows the statement for the system of Ito equations (f,g)- 
To prove the result for Stratonovich's interpretation we use the explicit relation between both 
interpretations and the representation ^ for the stochastic perturbations in the final part 
of the proof of Theorem pi This leads to the modified system of Ito equations (f,g), for 
which necessary and sufficient conditions are known. We observe that the conditions on the 
functions / and g are equivalent to the same conditions for the functions / and g, and are 
therefore invariant under the transformation relating both interpretations. □ 

Proof of Theorem 4- The comparison theorem for Ito's interpretation is valid by Theorem |8j 
To show the result for Stratonovich's interpretation of stochastic differential equations we 
use the explicit transformation formula, which leads to the modified system of Ito equations 
(f,g)- We apply the known result for Ito's interpretation and observe that the conditions for 
the functions / and g are equivalent to the conditions for the functions / and g. □ 

Concluding Remarks 

We obtained necessary and sufficient conditions for the invariance of rectangular subsets of 
the euclidean space under systems of stochastic differential equations and proved that the 
invariance property is independent of Ito's and Stratonovich's interpretation. In particular, 
we were able to characterize the class of stochastic perturbations that preserve the invariance 
property of the unperturbed deterministic system of ODEs. Such results are very relevant for 
applications and allow to validate stochastic models. 

When not only temporal but also spatial properties are relevant, the models are generally 
formulated as systems of stochastic PDEs. We are currently working on the extension of our 
invariance results for systems of parabolic PDEs under stochastic perturbations. A first result 
in this direction has been obtained in [1]. 
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